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Abstract 

FORM, a symbolic manipulation system, has been widely used in a lot of calculations for High 
Energy Physics due to its high performance and efficient design. Mathematica, another compu- 
tational software program, has also widely been used, but more for reasons of generality and 
user-friendliness than for speed. Especially calculations involving tensors and noncommutative 
operations like calculating Dirac traces can be rather slow in Mathematica, compared to FORM. 

In this article we describe FormLink and FeynCalcFormLink, two Mathematica packages to 
link Mathematica and FeynCalc with FORM. FormLink can be used without FeynCalc and 
FeynCalcFormLink, which is an extension loading FormLink and FeynCalc automatically. 

With these two packages the impressive speed and other special features of FORM get embedded 
into the generality of Mathematica and FeynCalc in a simple manner. 

FeynCalcFormLink provides a FORM-based turbo for FeynCalc, making it much more 
efficient. FormLink turns Mathematica into an editor and code organizer for FORM. 
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PROGRAM SUMMARY 

Title of program: FormLink/FeynCalcFormLink 



Available from: http://www.feyncalc.org/formlink/ 



Programming language: Mathematica 7, 8 or 9; C; FORM 4 
Computer: Any computer running Mathematica and FORM. 
Operating system: Linux, Windows, Mac OS X. 
External routines: FORM, C, FeynCalc. 

Keywords: Mathematica, MathLink, FORM, FeynCalc, ForniLink, FeynCalcForm- 
Link. 

Nature of physical problem: The functionality of FORM is restricted compared with 
Mathematica, which is broader, while its speed is slower for larger calculations involving 
tensors and noncommutative expansions. So how can we combine the speed of FORM 
with the versatility and broadness of Mathematical 

Method of solution: Named pipes, MathLink, C, Mathematica pre- and postprocessing. 

Typical running time: Seconds, minutes, or more, depends on the complexity of the 
calculation. There is some overhead for conversion of larger expressions. 
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LONG WRITE-UP 



I. Basic Ideas 

The basic and simplest way to communicate between Mathematica and FORM jjlfiGj is 
using input and output files. This method has been used in FeynCalc [7j and FormCalc 
[8], which prepares the symbolic expressions of the diagrams in an input file for FORM, 
runs FORM, and retrieves the results back to Mathematica. 

Another method to exchange data between different processes is to use pipes. The detailed 
usage of pipe communication between FORM and other external programs is described in 
[9]. In this article we implement a user- friendly communication between Mathematica and 
FORM using both possibilities. We would like to note that FormC ALC and FormGet . tm, 
mentioned at http://www.feynarts.de/formcalc, also use MathLink and pipes, but do 
not provide a general interface from Mathematica to FORM. 

The remainder of this section is intended for programmers only and not needed to un- 
derstand how to use FormLink. 

The basic idea is that we create two unnamed pipes, one is the read-only descriptor with 
file handler r# , the other is the write-only descriptor with file handler w#, then start FORM 
with the pipe option automatically through the MathLink executable FormLink 

form -pipe r#,w# init 

where init refers to the FORM code init.frm , which is discussed in the following. When 
the pipe connection has been established successfully, FORM sends its Process Identifer 
(PID) in ASCII decimal format with an appended newline character to the descriptor w# 
and then FORM will wait for the answer from the descriptor r#. The answer must be 
two comma-separated integers in ASCII decimal format followed by a newline character. 
The first integer corresponds to the FORM PID while the second one is the PID of parent 
process which started FORM. If the answer is not obtained after some time-out, or if it is 
not correct, i.e. it is not a list of two integers or the first integer is not the FORM PID, 
then FORM fails. When the channel has been established successfully, FORM will run the 
code in the init.frm file, containing the following instructions: 

Off Statistics; 
#ifndef 'PIPES.' 

#message "No pipes found"; 

. end; 
#endif 

#if ('PIPES.' <= 0) 

#message "No pipes found"; 

. end; 
#endif 

#procedure put(fmt, mexp) 

#toexternal 'fmt', 'mexp' 

#toexternal "#THE-END-MARK#" 
#endprocedure 
#setexternal 'PIPEl.'; 
#toexternal "OK" 
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#f romexternal 
. end 

The core parts are the last two hnes before the . end instruction. The first hne, 
#toexternal "OK", sends the word OK in ASCII string format from FORM to FormLink, 
it confirms that the communication channel has been established successfully, otherwise 
FormLink will treat it as failed. The second line blocks FORM and waits for the code 
which will be sent from Mathematica. When the prompt^ arrives, FORM will continue to 
execute code from #f romexternal. We defined a procedure named put to send data from 
FORM back to Mathematica. Note that if you want to send data without this procedure, 
you need to send the end mark, i.e. the string #THE-END-MARK#, to indicate that the data 
is complete, otherwise FormLink will be blocked until the end mark has been received. 

Until now we have demonstrated a round communication from Mathematica to FORM, 
and then back to Mathematica again. This procedure can be looped if we put another 
#f romexternal in the code sent from Mathematica to FORM, and in this sense, we get an 
interactive FORM, which cannot be easily achieved by the first way which exchanges data 
by input and output files. 



II. Installation 

To install the FormLink and FeynCalcFormLink packages, run the following in- 
struction in a Kernel or Notebook session of Mathematica 7, 8 or 9 on Linux, MacOSX, or 
Windows. 

Import ["http : //www. f eyncalc . org/f ormlink/ install .m"] 

The installer will automatically download formlink.zip from the urP and extract the files 
from the archive to the directory Applications in the directory $UserBaseDirectory, 
which is by default located in the search path of Mathematica. Their values for different 
platform are listed in TABLE |Tj Specifying 

$installdirectory = mydir 

before running the installer will change the installation directory. It is recommended to use 
a directory which is on the Mathematica path, e.g., HomeDirectory [] , or $BaseDirectory. 

For user convenience, the binary files of FormLink, form and tf orm for Linux, Microsoft 
Windows and Mac OS X are also installed. Furthermore the latest version of FeynCalc 



is downloaded from http : //www . f eyncalc . org| and installed automatically into the same 



directory unless it has been already installed somewhere on the Mathematica path. 
Platform $UserBaseDirectory 

Windows C : \Users\username\AppData\Roaming\Matheniatica 

Linux "/.Mathematica 
Mac OS X "/Library/Mathematica 

TABLE L The values of $UserBaseDirectory for different operating systems 



^The default prompt defined in FORM is a blank line, for details please see [S] 
^ttp : //www . f eynca lc . org/f ormlink/f ormlink.zip 



At the end of the installation FormLink and FeynCalcFormLink are loaded and two simple 
examples are run, one uses FormLink: 

FormLink [ " AutoDeclare vector p; 

Local T = g_(0, pl,p2,p3,p4,p5,p6) ; 
trace4 0; 
" (*, Form2M -> Identity *) 
(* <-- uncommenting returns a string *) 
IN- 
FORM and FormRead finished, time needed before translating to Mathematica: 0. sec 
Translation done. Total wall clock time needed: 0.1925 sec 
4 p 1 .p6 p2.p5 p3 .p4 - 4 p 1 .p5 p2.p6 p3 .p4 + 4 p 1 .p2 p3 .p4 p5 .p6 - 
4 p 1 .p6 p2.p4 p3 .p5 + 4 p 1 .p4 p2.p6 p3 .p5 + 4 p 1 .p5 p2.p4 p3 .p6 - 4 p 1 .p4 p2.p5 p3 .p6 + 
4 p 1 .p6 p2.p3 p4.p5 - 4 p 1 .p3 p2.p6 p4.p5 + 4 p 1 .p2 p3 .p6 p4.p5 - 4 p 1 .p5 p2.p3 p4.p6 + 
4 p 1 .p3 p2.p5 p4.p6 - 4 p 1 .p2 p3 .p5 p4.p6 + 4 p 1 .p4 p2.p3 p5 .p6 - 4 p 1 .p3 p2.p4 p5 .p6 

The other example executes FeynCalcFormLink, which generates and runs the corresponding 
FORM program, substituting ASCII values for greek indices intermediately: 

FeynCalcFormLink[ DiracTrace[GA[/./, v, p, a, z, a]]] 



AutoDeclare Index lor; 
Format Mathematica; 

L resFL = (g_(l,lor2)*g_(l,lor3)*g_(l,lor4)*g_(l,lor5)*g_(l,lor6)*g_(l,lorl)); 

trace4,l; 

contract 0; 

.sort; 

#call put("%E", resFL) 
#fromexternal 

Piping the script to FORM and running FORM 

Time needed by FORM : 0. seconds. FORM finished. Got the result back to Mathematica as a string. 
Start translation to Mathematica / FeynCalc syntax 

Total wall clock time used: 0.12 seconds. Translation to Mathematica and FeynCalc finished. 

4 glJO- gVp _4 ga(T gllT gVp _,_4 gOll g<TT gVp _4 gUT ^fl f gV <T _,_ 4 p gllT IT _,_ 

4 gl'P gVT -4 gCP gt"r gVT _,_4 gUT gpa- _^ ^av g/ir ^pa ^ p ^vt ^pcr _ 
4g"'^ gl'y gPT +4g"^' gt^o" gPT -^^gOifi gva- gPT _,_ 4 gap gpv gO-T _ ^ gav giip gO-T 

The source code of FormLink is included in the src folder. The scripts for compilation 
are listed in TABLE [TTj We used the GNU compiler gcc, producing an executable FormLink 
in the corresponding subdirectory in bin. 
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Platform 


make commands 




Microsoft Windows ("32 fe 64 bit) 


make 
make 


-f makefile, 
-f makefile, 


. cygwin 

. cygwin install 


Linux (32 bit) 


make 
make 


-f makefile, 
-f makefile. 


. Iinux32 

. Iinux32 install 


Linux (64 bit) 


make 
make 


-f makefile, 
-f makefile. 


. Iinux64 

. Iinux64 install 


Mac OS X (64 bit) 


make 
make 


-f makefile, 
-f makefile. 


.macosx64 
.macosx64 install 



TABLE II: The make commands for different operating systems. If needed, please modify the 
corresponding makefile, for example, change the location of the MathLink developer kit (MLINKDIR) 
if you did not use the default location during installation of Mathematica or if you want to use a 
different version than Mathematica 8. 



III. FormLink and FeynCalcFormLink Functions Reference 
A. Bgisic Usage 

The basic functions are FormLink to execute FORM programs written as a string in Math- 
ematica and FeynCalcFormLink to run a program specified in FeynCalc syntax through 
FORM. These two MathLink-hSiSed functions are easy to use for FORM and FeynCalc 
practicioners. 

• FormLink 

FormLink [" form statements "] runs the form statements in FORM and returns 
the result to Mathematica. If only one Local assignment is present the return value 
is the result. If more FORM Local variables are present, a list of results is returned. 
FormLink has several options, the default values follows, 

{Assign -> False, Form2M -> Form2M, FormSetup $FormSetup, Replace $Form2M, 
Print ^ True, Style ^ (RGBColor[0.444444, 0.222222, O.J, FontFamily ^ Courier}} 

Here we only explain the most important ones. 

The option Form2M can be used to specify the function which translates the ex- 
pression from FORM to Mathematica format, its default value is function Form2M, 
which has the following usage: Form2M [string, replist] translates string by 
ToExpression [StringReplace [string, replist] ,TraditionalForm] to Mathemat- 
ica. Form2M can also be set to Identity, which means that the result from FORM 
will be returned as a string, containing the exact output for the Local expression. If 
conversion to Mathematica is straighforward, i.e., there are no bracketed expressions 
or other syntaxes not easily interpreted by Mathematica, then it is best to use the set- 
ting Form2M— T-ToExpression, which will just call ToExpression on the string result 
coming from FORM. 

The option Replace is a list of string replacements, its default setting is $Form2M, 
containing a list of user-changeable (e.g. in Config.m) basic function-name translations 
from FORM to Mathematica. 
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{bernoulli_ -> BernouUiB, binom_ -> Binomial, cos_ Cos, fac_ -> Factorial, 
gcd_ GCD, ln_ -> Log, max_ -> Max, min_ Min, mod_ -> Mod, pi_ Pi, 
li2_ PolyLog2,, sqrt_ Sqrt, sign_ Sign, sin_ -> Sin, sum_ Sum, tan_ Tan, i_ -> I, 
e_ -> I*$LeviCivitaSign*Eps, d_ Pair, \ ^ , gi_ -> DiracGamma, ^ , _ ^ , [ ^ Hold[Identity][} 

The option Print can be used to switch informative messages on or off during exe- 
cution. The option FormSetup can be a hst of FORM settings, hke TempDir which 
dynamically produces a form. set file next to the binary which is then called auto- 
matically upon starting FORM. 

• FeynCalcFormLink 

FeynCalcFormLink [expr] translates the FeynCalc expression expr to a FORM 

program, identifying automatically traces, symbols, vectors, indices and the dimension, 
calculates it, pipes it back to Mathematica and translates it to FeynCalc syntax. 

There are several options for FeynCalcFormLink, the default settings are as follows, 
{Functions -> CFunctions, FeynCalcExtemal True, FormSetup :^ $FormSetup, 
Form2FC Form2FC, ExtraDeclare {}, IDStatements -» {}, Print True, 
Replace ^ {}, Style ^ {RGBColor[0.444444, 0.222222, 0.], FontFamily ^ Courier} 

Most of the options are similar to FormLink, we only explain Functions, 
ExtraDeclare, IDStatements and Form2FC. 

If the option Functions is set to "CFunctions", then all non-System' functions, 
except those present in $M2Form and some FeynCalc functions, are automatically 
declared CFunctions in FORM. If Functions option is set to "Functions", then they 
are declared noncommutative functions, i.e.. Functions in FORM. 

The option ExtraDeclare can be used to put any extra valid FORM declarations 
which are not identified automatically, for example: 

ExtraDeclare— 7'{" CFunctions GammaFunction; " , "Functions MyOperator ; "} 

The option IDStatements can be set to a string or a list of strings corresponding to 
FORM identify statement like {"id kl .kl=mass"2; "}. 

The option Form2FC is set to the function being used to translate expression from 
FORM to FeynCalc format, the default is Form2FC. 



B. Advanced Usage 

The basic internal procedure to use FormLink is to start FORM with FormStart, 
then send the code to FORM for execution using FormWrite and FormPrompt, read and 
convert the result back to Mathematica through FormRead, finally, stop FORM by caUing 
FormStop. All these calls can be encoded into a single function named FormLink which has 
been introduced in the previous section. 

• FormStart 

FormStart [] automatically determines the path of the form executable, which is lo- 
cated in the corresponding subfolder in the bin directory, and then starts FORM in 
pipe mode. You can also call FormStart [f ormpath] with explicit full path f ormpath 
of the form excutable. 
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• FormWrite 

FormWrite [script] sends the script with an appended newhne character, to 
FORM. Notice that FORM will not start to execute the script right away and you 
can send your scripts by using the FormWrite [script] many times. When you are 
finished, send the prompt, then FORM executes all instructions sent. 

FormWrite [scripts] , where scripts is a list of strings, will call FormWrite [script] 

for each element script in the list scripts, so if your code spans several lines, you can 
also put them as a list with each element corresponding to a single line of your scripts. 

• FormPrompt 

FormPrompt [] sends the prompt to FORM, which will make FORM continue to exe- 
cute the code you have sent. The default prompt in FORM is a blank line, FormLink 
has adopted this default option, so do not send blank lines unintentionally. 

• FormRead 

FormRead[] reads the data from FORM. It should be noted that if no data is sent 
from FORM, the calling thread will be blocked until data is available. You can use 
the procedure put defined in init .f rm to send the data from FORM. 

FormRead [] will first check whether the pipe has been closed or not, if the pipe has 
been closed, for example when FORM has encountered some problems, FormRead [] 
will redirect the standard output from FORM to Mathematica, so the user can check 
the error messages. 

• FormStop 

FormStop[] uninstalls the link to Form. FormStop [All] kills all running FormLink 
processes. 

As we discussed in section |T| there are two ways to communicate between FORM and 
Mathematica, the functions introduced above are all used with the method of piping. We 
also provide two functions which are not using MathLink, but deal with input and output 
files only. 



• RunForm 

RunForm [script] runs script in FORM and writes the result to runform.frm and 
the log file to form. log in the current directory, i.e., the value returned by the 
Mathematica function Directory [] . 

RunForm [script , f ormf ile] uses f ormf ile instead of runform.frm. The first argu- 
ment script can be a string or a list of strings. 

An optional third argument can be given to use a specific FORM executable, otherwise 
a FORM executable from $FormLinkDir/bin is used. 

• ReadString 

ReadString[str] imports str as Text and translates it to Mathematica syntax by 
using $Form2M. You can use #write preprocessor to write your data to the output file 
by FORM, and use ReadString to read it into Mathematica. 
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To facilitate the usage of FormLink with FeynCalc, two functions are provided for 
performing the conversions between FeynCalc and FORM^: 

• FC2Form 

FC2Form[exp] translates exp in FeynCalc format to FORM format. FC2Form [exp] 

returns a list of two elements, the first one is a list containing the script which will 
be sent to FORM. You can use the function ShowScript [script] to display the 
script. The second one is a list of replacement rules, which will be used in Form2FC 
to translate the result from FORM format back to FeynCalc. FC2Form has the 
options : Functions, Dimension, ExtraDeclare, IDStatements, Print and Replace: 

{Functions -> CFunctions, Dimension -> Automatic, 
ExtraDeclare {}, IDStatements -> {}, Print -> True, Replace -> {} 

• Form2FC 

Form2FC [f ormexpr] is used to translate f ormexpr which is in FORM format back to 
the FeynCalc one. Form2FC[f ormexpr, replacelist] apphes the substitution hst 
replacelist at the end. The second argument is usually the second item of of the 
list returned by FC2Form. 

So the general steps to use the package with FeynCalc are to convert FeynCalc code 

to FORM, and then send the converted code to FORM for executing with FormLink, 
and finally convert the results in FORM format back to FeynCalc. All these steps are 
implemented into a single function FeynCalcFormLink, which has been introduced in the 
previoius section. 

IV. Examples Using FormLink, FeynCalcFormLink and RunForm 

We list a few examples using the FormLink, FeynCalcFormLink and RunForm functions. 
More examples can be found in the Examples directories of both packages. 

A. Using FormLink 

FormLink is a function for running FORM from Mathematica and returning the result 
to Mathematica. 

A short trace 
short = FormLink [ " 

AutoDeclare vector p; 

Local T = g_(0, pl,p2,p3,p4,p5,p6) ; 

trace4 ; 

"] 



^There are also two older FeynCalc functions named FeynCalc2F0RM and F0RM2FeynCalc to perform 
the conversions. 
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FORM and FormRead finished, time 

needed before translating to Mathematica: 0.001 sec 

Translation done. Total wall clock time needed: 0.124 sec 

4 pl.p6 p2 .p5 p3.p4 - 4 pl.p5 p2.p6 p3.p4 - 4 pl.p6 p2.p4 p3.p5 + 
4 pl.p4 p2 .p6 p3.p5 + 4 pl.p5 p2.p4 p3.p6 - 4 pl.p4 p2.p5 p3.p6 + 
4pl.p5p2.p3p4.p5-4pl.p3p2.p6p4.p5 + 4pl.p2p3.p6p4.p5- 
4pl.p5p2.p3p4.p5 + 4pl.p3p2.p5p4.p6-4pl.p2p3.p5p4.p6 + 
4 pi .p4 p2 .p3 p5.p5 - 4 pl.p3 p2.p4 p5.p6 + 4 pl.p2 p3.p4 p5.p6 

% 1 1 InputForm 

4 * pi . p6 * p2 . p5 * p3 . p4 - 



4 * pi 


. p5 * p2 


. p6 * p3 


. p4 - ^ 


\ * pi 


. p6 * p2 


. p4 ■ 


* p3 


. p5 + 


4 * pi . 


p4 * p2 . 


p6 * p3 . 


p5 + 4 


* pi . 


p5 * 


p2 . 


p4 * 


p3 . 


p6 - 


4 * pi 


. p4 * p2 


. p5 * p3 


. p6 + 














4 * pi . 


p6 * p2 . 


p3 * p4 . 


p5 - 4 


* pi . 


p3 * 


p2 . 


p6 * 


p4 . 


p5 + 


4 * pi 


. p2 * p3 


. p6 * p4 


. p5 - 














4 * pi . 


p5 * p2 . 


p3 * p4 . 


p6 + 4 


* pi . 


p3 * 


p2 . 


p5 * 


p4 . 


p6 - 


4 * pi 


. p2 * p3 


. p5 * p4 


. p6 + 














4 * pi . 


p4 * p2 . 


p3 * p5 . 


p6 - 4 


* pi . 


p3 * 


p2 . 


p4 * 


p5 . 


p6 + 


4 * pi 


. p2 * p3 


. p4 * p5 


. p6 















Using the option Form2M— j-Identity returns a string: 
short = FormLink[ " 

AutoDeclare vector p; 
Local T = g_(0, pl,p2,p3,p4,p5,p6) ; 
trace4 0; 
", Form2M -» Identity 

] 

FORM and FormRead finished, time 

needed before translating to Mathematica: 0.001 sec 
Translation done. Total wall clock time needed: 0.1185 sec 

(4*pl .p2*p3 .p4*p5 .p6-4*pl .p2*p3 .p5*p4 .p6+4*pl .p2*p3 .p6*p4 .p5-4*pl .p3*p2 .p4 

*p5 .p5+4*pl .p3*p2 .p5*p4 .p5-4*pl .p3*p2 .p6*p4 .p5+4*pl .p4*p2 .p3*p5 .p5-4*pl . 

p4*p2 .p5*p3 .p5+4*pl .p4*p2 .p6*p3 .p5-4*pl .p5*p2 .p3*p4 .p5+4*pl .p5*p2 .p4*p3 . 

p5-4*pl .p5*p2 .p5*p3 .p4+4*pl .p5*p2 .p3*p4 .p5-4*pl .p5*p2 .p4*p3 .p5+4*pl .p5* 

p2 .p5*p3 .p4 ) 
short // StringQ 
True 

Clear [short] 
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A longer trace 
meditim = FormLink [ " 

AutoDeclare vector p; 
Local T = g_(0, 
Pl,p2,p3,p4,p5,p6,p7,p8,p9,pl0,pll,pl2,pl3,p3) ; 
trace4 0; 

II 

r 

Form2M :^ Function[x, ToExpression [StringReplace [x, " . " -» " 

(* Since we want to change pl.p2 to plp2, 
we replace the . to an empty string "*) 

]; 

FORM and FormRead finished, time 

needed before translating to Mathematica: 0.0595 sec 
Translation done. Total wall clock time needed: 0.3055 sec 
mediiun // Length 

11 763 

mediTom [ [ 1 ] ] 

4 pl0pl3 pllpl2 plp9 p2p8 p3p3 p4p7 p5p6 
Variables [medium] // InputForm //Short 

{pl0pl3, pllpl2, plp9, p2p8, p3p3, p4p7, 

p5p6, plOpl2, << 67 >>, p9pl0, p9pll, p9pl2, p9pl3} 
Clear [medivim] ; 

Special cases 

Objects like [1-x] in FORM are translated to (1-x) in Mathematica 
FormLink [" 

Symbols y, [1-x] ; 
L F = [1-x] * (2*y-l) '^2*3; 

"] 

FORM and FormRead finished, time 

needed before translating to Mathematica: 
Translation done. Total wall clock time needed 
3 (1-x) -12 (1-x) y + 12 (1-x) 
% / / Factor 
-3 (-1 +x) (-1 +2 y)2 
FormLink [I] 

FORM and FormRead finished, time 

needed before translating to Mathematica: 
Translation done. Total wall clock time needed 
i 



.0005 sec 
0.1075 sec 



,001 sec 
0.0855 sec 
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Polynomial examples 
Expand {a + b)^ 

FormLink[" Symbols a,b; Local F = (a+b)^2"] 

FORM and FormRead finished, time 

needed before translating to Mathematica: 0.0005 sec 
Translation done. Total wall clock time needed: 0.085 sec 

a^ + 2 a b + 

Alternatively we may just enter the Mathematica expression which is then automatically 
translated to a FORM program: 

FormLink[ (a + b) '^2] 

FORM and FormRead finished, time 

needed before translating to Mathematica: 0.001 sec 
Translation done. Total wall clock time needed: 0.0855 sec 

a^ + 2 a b + b^ 

Expand {a + b + c + d + e + f + gf'^ 

The default option setting of Form2M does general translation of FORM syntax to Mathe- 
matica /FevnCalc. Here we can use the simpler ToExpression as a setting of Form2M. 
While it is possible to use even 42 instead of 21, we see already for power 21 that it takes 
longer to transfer the large expression between FORM and Mathematica than to calculate 
it in either FORM or Mathematica directly. So it is best not to transfer large expressions, 
if possible. 

fresult = 
FormLink [ " 

Symbols a, b, c, d, e, f , g; 

Local F = (a+b+c+d+e+f +g) '^21; 

n 

r 

Form2M -» ToExpression 

]; 

FORM and FormRead finished, time 

needed before translating to Mathematica: 1.6787 sec 
Translation done. Total wall clock time needed: 5.9193 sec 
fresult // Length 

296 010 

So for somewhat larger expressions there is noticable overhead for piping the expression 
back to Mathematica as a string by FormRead and for translating the string to Mathematica 
syntax (ToExpression). 

An accurate timing of FORM can be done by using RunForm: 
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RunForm®" On ShortStats; 

Symbols a, b, c, d, e, f , g; 

Local F = (a+b+c+d+e+f +g) '^21; 
II . 

r 

FORM 4 . (Aug 31x2012) 32 - bits Run : Wed Dec 12x12 : 01 : 46x2012 



On ShortStats; 




Symbols a. 


b, c, d, e, f, g; 




Local F = 


(a + b + c + d + e + f + g) '"21; 




. end 






0.24s 


1 > 100 000 100 000 : 


3 015 192 F 


.53 s 


1 > 200 000 200 000 : 


6 000 862 F 


0.81s 


1> 296 010 296 010: 


8 643 114 F 


0.85s 


1 > 296 010 -- > 296 010 


: 8 643 07 


0.85 sec out of 0.85 sec 





8 F 



Expand {a + h + c + d + e + f + gY'^ using RunForm 

Last [StringSplit [7^ ,"\n"] ] & is a pure function to extract the last line of the string 
output of RunForm, returning the timing. 

Last [StringSplit[tt, "\n"]] &@ 
RunForm[" On ShortStats; 

Symbols a, b, c, d, e, f , g; 
Local F = (a+b+c+d+e+f+g) *42; 
" , Print -* False] 
83.85 sec out of 86.46 sec 
First @ AbsoluteTiming[ mresult = Expand[ (a + b + c + d + e + f + g) ^ 42] ; ] 
63.248618 

{Length® mresult, MemorylnUse [ ] , Clear [mresult] , MemorylnUse [ ] } 

{ 1 2 27 1 5 1 2, 2 039 8 1 7 984, Null, 1 65 005 1 84} 

So mresult has more than 12 million terms and it occupies around 2 GB. For this example 
with a lot of terms FORM 4 is slightly slower than Mathematica 9. 

Clear [f result, mresult] 

(a + h)^ and a% — )• c 
FormLink [ " 

Symbols a,b,c; 
Local F=(a+b)'^6; 
id a''2*b = c; 

"] 

FORM and FormRead finished, time 

needed before translating to Mathematica: 0.0005 sec 
Translation done. Total wall clock time needed: 0.217 sec 

a'^ + 6ab^+b'^ + 6a^c + 20ab2c + 15b^c+15c^ 



13 



To do the same in Mathematica requires more work: 

Expand[ (a + b) ^] // . 
(a'' (n_Integer /; n > 1) * (b-" (m_. ) ) ^ ( (a'^ (n - 2) c) b'' (m-1))) 

a^ + 6ab^+b^ + 6a^c + 20ab^c + 15b^c + 15c^ 

Multiple Local variables: Local Fl= {a + b + c)^° Local F2 = (a + 6 + c + 
Length /@ FormLink [ " 

tt:SmallSize 2000; 
tt:LargePatches 4; 
Symbols a,b,c,d; 
Local Fl = (a+b+c) '^10; 
Local F2 = (a+b+c+d) ''lO; 

"] 

FORM and FormRead finished, time 

needed before translating to Mathematica: 0.002 sec 
Translation done. Total wall clock time needed: 0.1375 sec 
{66, 286} 

Length/® Expand [{ (a + b + c)^°, (a + b + c + d)"}] 

{66, 286} 

Commuting and noncommuting functions 

FormLink [ " 

Functions A1,B1; 

CFunctions A2 , B2 ; 

Local Fl = (A1+B1)'"3; 

Local F2 = (A2+B2)'^3; 
", Form2M M Function [x, ToExpression[StringReplace[x, "*"^ "**"]]] 
] 

FORM and FormRead finished, time 

needed before translating to Mathematica: 0.001 sec 
Translation done. Total wall clock time needed: 0.0885 sec 
{a1 **A1 **A1 +A1 **A1 **B1 +A1 **B1 **A1 +A1 **B1 **B1 + B1 **A1 **A1 + 

B1**A1**B1+B1**B1**A1+B1**B1**B1, A2^+B2-^ + 3**A2**B2^ + 3**A2^**B2| 

Index and Vector (Local F=pl(il)*(p2(il)+p3(i3))*(pl(i2)+p2(i3));) 
tmp = FormLink [ " 

Index il,i2,i3; 
Vector pl,p2,p3; 

Local F=pl (il) * (p2 (il) +p3 (13) ) * (pi (12) +p2 (13) ) ; 

"] 

FORM and FormRead finished, time 

needed before translating to Mathematica: 0.001 sec 
Translation done. Total wall clock time needed: 0.093 sec 
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p2.p3pl[il] +pl.p2pl[i2] +pl.p2p2[i3] +pl[il] pl[i2] p3[i3] 
tmp // InputForm 

p2 . p3*pl[il] + pi . p2*pl[i2] + pi . p2*p2[i3] + pl[il] *pl[i2] *p3[i3] 
tmp = FormLink [ " 

Index 11,12,13; 
Vector pl,p2,p3; 

Local F=pl (11) * (p2 (11) +p3 (13) ) * (pi (12) +p2 (13) ) ; 

" , Form2M ^ Identity] 

FORM and FormRead finished, time 

needed before translating to Mathematica: 0.0005 sec 
Translation done. Total wall clock time needed: 0.091 sec 
(pi (il) *pl (12) *p3 (13) +pl (il) *p2 .p3+pl (12) *pl .p2+p2 (13) *pl .p2) 
tmp // InputForm 

(pi (il) *pl (12) *p3 (13) +pl (11) *p2 .p3+pl (12) *pl . p2+p2 ( 13 ) *pl . p2 ) 
Needs [ "FeynCalcFormLlnk" " ] 
tmp2 = 

Form2FC[tmp, {pi Momentiim[pl] , p2 Momentum [p2 ] , p3 -» Momentum[p3] , 
11 -» LorentzIndex[ll] , 
12 -» Lorentz Index [12] , 13 -» LorentzIndex[13] } ] 

pi" pV^ p3" + pi" p2 • p3 + pl'^ pi • p2 + p2" pi • p2 
tmp2 // InputForm 

FV[pl, il]*FV[pl, i2]*FV[p3, 13] + 
FV[pl, i2]*SP[pl, p2] + FV[p2, i3]*SP[pl, p2] + 
FV[pl, 11] * SP [p2, p3] 
tmp2 // FCI // InputForm 

Pair [Lorentzlndex [11 ] , Momentum [pi ] ] * 
Pair [Lorentzlndex [12 ] , Momentum [pi ] ] * 
Pair [Lorentzlndex [ 13 ] , Momentum [p3 ] ] + 
Pair [Lorentzlndex [12 ] , Momentum [pi ] ] * 

Pair [Momentum [pi ] , Momentum [p2 ] ] + Pair [Lorentzlndex [ 13 ] , Momentum [p2 ] ] * 
Pair [Momentum [pi ] , Momentum [p2 ] ] + Pair [Lorentzlndex [ 11 ] , Momentum [pi ] ] * 
Pair [Momentum [p2 ] , Momentum [p3 ] ] 
Clear [tmp, tmp2] 

This sets back the default output format type to StandardForm. 
SetOptlons [tt, "CommonDef aultFormatTypes" -» 

{"Input" StandardForm, " Inputlnllne" -» StandardForm, 
"Output" -* StandardForm, "Outputlnllne" -» StandardForm, 
"Text" -» TextForm, "Textlnllne" -» TradltlonalForm} ] & /@ 
{$FrontEnd, $FrontEndSesslon} ; 
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B. Using FeynCalcFormLiiik 



Examples 1 

Calculate a trace tr(7''7'^7''7'^7'^7°"). Using DiracTrace in FeynCalc does not calculate 
immediately: 

spur = DiracTrace [GAD [/u, v, p, a, z, a]] 

Feeding this into FeynCalcFormLink has enough information to tell FORM to do the trace 
in D dimensions: 

FeynCalcFormLink [spur] 



Symbol D; 
Dimension D; 
AutoDeclare Index lor; 
Format Mathematical 

L resFL = (g_(l,lorl)*g_(l,lor2)*g_(l,lor3)*g_(l,lor4)*g_(l,lor5)*g_(l,lor4)); 

tracen,!; 

contract 0; 

.sort; 

#call put("%E", resFL) 
#fromexternal 

Piping the script to FORM and running FORM 

Time needed by FORM : 0.006 seconds. FORM finished. Got the result back to Mathematica as a string. 
Start translation to Mathematica / FeynCalc syntax 

Total wall clock time used: 0.15 seconds. Translation to Mathematica and FeynCalc finished. 

- 4 D ^ P + 4 D g*^ P g"^ " - 4 D gf ^ " + 8 " P - 8 g'^ P g"^ ^ + 8 ^ g'' ^ 
Factor@% 

-4 (D - 2) (g'^-g'P - g^Pg'- + gP'gn 

FeynCalcFonnLink[DiracTrace[GA[a, ^, a, z, 5]]] 



AutoDeclare Index lor; 
Format Mathematica; 

L resFL = (g_(l,lorl)*g_(l,lor2)*g_(l,lor3)*g_(l,lor4)*g5_(l)); 

trace4,l; 

contract 0; 

.sort; 

#call put("%E", resFL) 
#fromexternal 

Piping the script to FORM and running FORM 

Time needed by FORM : 0.003 seconds. FORM finished. Got the result back to Mathematica as a string. 
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Start translation to Mathematica / FeynCalc syntax 

Total wall clock time used: 0.12 seconds. Translation to Mathematica and FeynCalc finished. 

^ gO /3 a t 

Notice that loading FeynCalcFormLink by default sets $LeviCivitaSign=-I, such that 
traces involving 7^ agree. 

$LeviCivitaSign 

-i 

FeynCalcFormLink [ 

fun[bla] DiracTrace [somefunction[/3] anotherf unction [0] GK{n, v, p, a, 7]], 
Print -» True] // Factor 



Symbols bla,syml,sym2; 

AutoDeclare Index lor; 

AutoDeclare Symbol sym; 

CFu notions anotherfunction,fun,somef unction; 

Format Mathematica; 

L resFL = (anotherfunction(sym2)*fun(bla)*g_(l,lorl)*g_(l,lor2)*g_(l,lor3)*g_(l,lor4) 
*(g7_(l)/2)*somefunction(syml)); 

trace4,l; 
contract 0; 
.sort; 

#call put("%E", resFL) 
#fromexternal 

Piping the script to FORM and running FORM 

Time needed by FORM : 0.003 seconds. FORM finished. Got the result back to Mathematica as a string. 
Start translation to Mathematica / FeynCalc syntax 

Total wall clock time used: 0.18 seconds. Translation to Mathematica and FeynCalc finished. 

- 2 anotherf unction (0) fun(bla) 

somefunction(/3) (- g^^ " g'^ ^ + g^"- ^ g"" - g^"" g'"' + e^"""") 
TR[GA[Ai, V, p, a, 7]] //Factor 

Example 2 

FeynCalcFormLink [DiracTrace [ (m + GS [p] ) . (M + GS [q] ) ] ] 



Symbols m,M; 

Vectors p,q; 

Format Mathematica; 

L resFL = ((m*gi_(l)+g_(l,p))*(M*gi_(l)+g_(l,q))); 

trace4,l; 

contract 0; 

.sort; 

#call put("%E", resFL) 
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#fromexternal 



Piping the script to FORM and running FORM 

Time needed by FORM : 0.002 seconds. FORM finished. Got the result back to Mathematica as a string. 
Start translation to Mathematica / FeynCalc syntax 

Total wall clock time used: 0.12 seconds. Translation to Mathematica and FeynCalc finished. 

4 m M + 4 p • q 

Example 3 

Define simple typesetting rules for pi, ... p8 : 

Do [With [{v = Symbol ["p" <>ToString[i] ] }, 

MakeBoxes[v, TraditionalForm] = SubscriptBox[ "p" , i] ] , {i, 8}] 

Enter a trace tr[(|^i + m)j^{^2 + |^ + "^)7^t|^2] like this in FeynCalC: 

exp = DiracTrace[ (GS[pl] + m) .GAD[^] . (GS [p2 + k] + m) .GhD[^] .GS[p2] ] 

tr((m + y ■ pl).y^.(y ■ {k + pj) + m).y^.{y ■ p2)) 

Calculate it through FORM: 
Rl = FeynCalcFonnLink@exp 



Symbols D,m; 
Dimension D; 
Vectors k,pl,p2; 
AutoDeclare Index lor; 
Format Mathematica; 

L resFL = ((m*gi_(l)+g_(l,pl))*g_(l,lorl)*(m*gi_(l)+g_(l,k)+g_(l,p2))*g_(l,lorl) 

*g_(i,p2)); 

tracen.l; 
contract 0; 
.sort; 

#call put("%E", resFL) 
#fromexternal 



Piping the script to FORM and running FORM 

Time needed by FORM : 0.002 seconds. FORM finished. Got the result back to Mathematica as a string. 
Start translation to Mathematica / FeynCalc syntax 

Total wall clock time used: 0.13 seconds. Translation to Mathematica and FeynCalc finished. 

-ADmk-p2-^Dmp2 + 'iDmpi ■ p2 + 8 m k ■ p2 + Q m p\ 

Check with direct calculation in FeynCalC: 
R2 = exp / . DiracTrace -» Tr 
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4i^-D mk ■ p2 - Dm p\+ Dm pi ■ p2 + 2mk ■ p2+2m pj) 

Expand® FCE [Rl - R2] 



Clear [Rl, R2] 

C. Using RunForm 

RunForm runs a FORM program without MathLink, but by calling FORM directly. 
Opt ions ©RunForm 

{Style ^ (RGBColor[0.444444, 0.222222, 0.], FontFamily ^ Courier}, 

FormSetup :-> $FormSetup, Print -> True} 
SetDirect ory [ $TemporaryDirectory] ; 

AbsoluteTiming [ 
RunForm [ " symbol a,b,c; 

Format NoSpaces; 
On ShortStats; 
L res= (a+b+c)''5; 
. sort 

ttwrite <res.txt> \"(%E)\", res 
.end 

"tmp.frm"] ; 
result = Form2M[ReadString["res.txt"] ] ; 
result] 



FORM 4.0 (Aug 31 2012) 32-bits 
symbol a,b,c; 
Format NoSpaces; 
On ShortStats; 
L res= (a+b+c)^5; 
.sort 

0.00s 1> 21-> 

#write <res.txt> "(%E)", res 

.end 

0.00s 21> 21-> 

0.00 sec out of 0.00 sec 



{0.147502, a^ + Sa^b+Sa" c + 10 + 20 b c + 

1 a^ + 10 a^ b^ + 30 a^ b^ c + 30 a^ b + 10 a^ + 5 a i?" + 2 a c + 
30 a b^ + 2 a b + 5 a c" + + 5 c + 10 b^ + 10 b^ + 5 b c" + } 

FilePrint ["res .txt "] 



Run: Wed Dec 12 13:21:06 2012 

21: 380 res 

21: 380 res 
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(c^5+5*b*c^4+10*b^2*c^3+10*b^3*c^2+ 

5*b^4*c + b'"5 + 5*a*c''4+20*a*b*c''3 + 30*a* 
b'~2*c''2+20*a*b'"3*c + 5*a*b^4 + 10*a^2*c''3 + 
30*a''2*b*c^2 + 30*a^2*b'^2*c + 
10*a''2*b''3 + 10*a''3*c-^2 + 20*a'^3*b*c + 
10*a'"3*b'^2 + 5*a'"4*c + 5*a''4*b + a'"5) 

SetDirectory [$TemporaryDirectory] ; Absolu'teTiining[ 



RunForm [ { 

" symbol a,b,c; 

Off Statistics; 
L res= (a+b+c) '"199; 
. sort; 

ttwrite <res.txt> \"(%E)\", res 

II 

}, "tmp.frm"]; 
result = Form2M[ReadString["res.txt"] ] ; 
result // Length] 



FORM 4.0 (Aug 31 2012) 32-bits 
symbol a,b,c; 
Off Statistics; 
L res= (a+b+c)^199; 
.sort; 

#write <res.txt> "(%E)", res 
.end 

1.21 sec out of 1.22 sec 



Run: Wed Dec 12 13:21:06 2012 



{5 . 324083, 20 100} 

Length [Expand [ (a + b + c) 199] - result] 



SetDirectory [$TemporaryDirectory] ; AbsoluteTiming[ 
RunForm [ { 

" symbol a,b,c; 

Off Statistics; 

L res= (a+b+c) '^199; 

. sort; 

ttwrite <res.txt> \"(%E)\", res 

}, "tmp.frm"]; 
(*result = Form2M[ReadString["res.txt"] ] ; *) 

(* if the result from FORM is easy to load into Mathematica, 
Get is faster : *) 
result = Get["res.txt"] ; 
result // Length] 



FORM 4.0 (Aug 31 2012) 32-bits Run: Wed Dec 12 13:21:12 2012 
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symbol a,b,c; 

Off Statistics; 

L res= (a+b+c)^199; 

.sort; 

#write <res.txt> "(%E)", res 
.end 

1.21 sec out of 1.21 sec 



{1 . 605025, 20 100} 

Length [Expand [ (a + b + c) ^199] - result] 



For such simple algebraic operations Mathematica is actually much faster. 
AbsoluteTiming [Expand [ (a + b + c) ''199] ; ] 

{0.072501, Null} 

SetDirectory [$TemporaryDirectory] ; 
AbsoluteTiming [ 
RunForm [ { 

"symbol a, b, c, d, e, f , g; 
Format NoSpaces; 
Off Statistics; 
L res= (a+b+c+d+e+f +g) "21; 
. sort 

ttwrite <res.txt> \"(%E)\", res" 
}, "tmp.frm"]; 

result = Get["res.txt"] ; (*Form2M [ReadString[ "res . txt" ] ] ; *) 
result // Length] 



FORM 4.0 (Aug 31 2012) 32-bits Run: Wed Dec 12 13:21:14 2012 

symbol a,b,c,d,e,f,g; 
Format NoSpaces; 
Off Statistics; 

L res= (a+b+c+d+e+f+g)^21; 
.sort 

T^write <res.txt> "(%E)", res 
.end 

1.34 sec out of 1.35 sec 

{5 . 655087, 296 010} 

AbsoluteTiming [mresult = Expand[ (a + b + c + d + e + f + g) ''21] ; ] 
result - mresult 

{1.003515, Null} 




21 



V. APPLICATION TO TREE LEVEL PROCESSES 



A. e+e — 7- r+T ud ^v^VrVr 

Let us consider some applications of FormLink combined with FeynCalc. We take the 
process: e^e~ — )■ r+r^ — t- ud ^v^Vt-Vt as an example, which has been considered in the 
Form courses[10j. We can express the squared amplitude as: 

^ Tr[(;^2 - me)7^HA + "^e)7"^] 

* Tr[(;^4 + m4)7''2 77(^5 - rn^)^"'^ ^^^^ 

* Tr[(;^7 + m7)7''3^7(|^8 - "1.8)7'''' 77] (1) 

The expression in ([T| can be easily translated to FeynCalc syntax and then executed by 
FeynCalcFormLink, it takes a fraction of a second to run the code. The program is: 

incomingemeppalr = DiracTrace[ (GS[p2] -emass) .GA[inl] . 

(GS [pi] + emass) .GA[nl] ] ; 
tauline = DiracTrace[ (GS[p3] +mass3) .GA[m2] .GA[7] . 

(GS[ql] +tmass) .GA[ml] . 

(-GS [q2] + tmass) . 

GA[m3] .GA[7] . (GS[p6] -mass6) . 

GA[n3] .GA[7] . (-GS[q2] + tmass) .GA[nl] . 

(GS[ql] +tmass) .GA[n2] .GA[7] 

]; 

udbarpair = DiracTrace[ 

(GS[p4] +mass4) .GA[m2] .GA[7] . (GS[p5] -mass5) .GA[n2] .GA[7] ] ; 
nubarmupair = 

DiracTrace[ (GS[p7] +mass7) .GA[m3] .GA[7] . (GS[p8] -mass8) .GA[n3] .GA[7] 

]; 

1 

formexample = Incomingemeppalr. tauline. udbarpair. nubarmupair 

2I6 

tr((-y • p2 - emass).7™'. (emass + y ■ pl).-y"'). 

256 ^ 

tr((mass3 + y ■ ^3>).y"''^ .y"^ .(y ■ ql + tmass). y™'. (tmass - y ■ q2).y^^. 

y^ -{y ■ p6 - mass6).7"^.'y''. (tmass - y ■ q2).7"'.(7 • ql + tmass). 7"^ .7^). 
tr((mass4 + 7 • p4).7™^.7''.(7 • p5 - mass5).7"^.7^).tr((mass7 + 7 • p7).7™^7''.(7 • p8 - mass8).7"^.7'') 
Options [FeynCalcFormLink] 

{Functions CFunctions, FeynCalcExternal True, FormSetup :^ $FormSetup, 
Form2FC Form2FC, ExtraDeclare {}, IDStatements -> {), Print -> True, 
Replace ^ (}, Style ^ {RGBColor[0.444444, 0.222222, 0.], FontFamUy ^ Courier)) 
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R3 = FeynCalcFonnLink[f ormexample, IDStatements -» { 
"id ql.ql = tmass'^2; 
id q2.q2 = tmass''2; 
id pl.p2 = s/2-emass'^2; 
id ql.q2 = s/2-tmass''2 ; 

n 

}] // Function[x, Collect2[x, {tmass, emass}]] 



Symbols emass, mass3,mass4,mass5,mass6,mass7,mass8,s, tmass; 

Indices ml,m2,m3,nl,n2,n3; 

Vectors pl,p2,p3,p4,p5,p6,p7,p8,ql,q2; 

Format Mathematical 

L resFL = (((-(emass*gi_(l))+g_(l,p2))*g_(l,ml)*(emass*gi_(l)+g_(l,pl)) 

*g_(l,nl)*(mass4*gi_(2)+g_(2,p4))*g_(2,m2)*(g7_(2)/2)*(-(mass5*gi_(2)) 

+g_(2,p5))*g_(2,n2)*(g7_(2)/2)*(mass7*gi_(3)+g_(3,p7))*g_(3,m3) 

*(g7_(3)/2)*(-(mass8*gi_(3))+g_(3,p8))*g_(3,n3)*(g7_(3)/2)*(mass3 

*gi_(4)+g_(4,p3))*g_(4,m2)*(g7_(4)/2)*(tmass*gi_(4)+g_(4,ql)) 

*g_(4,ml)*(tmass*gi_(4)-g_(4,q2))*g_(4,m3)*(g7_(4)/2)*(-(mass6 

*gi_(4))+g_(4,p6))*g_(4,n3)*(g7_(4)/2)*(tmass*gi_(4)-g_(4,q2)) 

*g_(4,nl)*(tmass*gi_(4)+g_(4,ql))*g_(4,n2)*(g7_(4)/2))/256); 

trace4,l; 

trace4,2; 

trace4,3; 

trace4,4; 

contract 0; 

.sort; 

id ql.ql = tmass^2; 
id q2.q2 = tmass^2; 
id pl.p2 = s/2-emass^2; 
id ql.q2 = s/2-tmass^2; 
.sort; 

#call put("%E", resFL) 
#fromexternal 

Piping the script to FORM and running FORM 

Time needed by FORM : 0.017 seconds. FORM finished. Got the result back to Mathematica as a string. 
Start translation to Mathematica / FeynCalc syntax 

Total wall clock time used: 0.2 seconds. Translation to Mathematica and FeynCalc finished. 
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4 emass^ tmass^ pa ■ Pi Pe • Pa 

{-2 p-j • ql p5 • q2 - 2 p^ ■ ql p-j • q2 - 2 p^ • ql pv • ql - 2 ps • q2 pv • q2 + s ps • Pv ) + 
8 emass^ s P3 ■ p^ p^ ■ pg ps • ql p7 • q2 - 

2 tmass^ P3 • P4 ps • Ps ( - 2 s p? • ql ps • q2 - 2 s ps • ql pv • q2 - 

4 ps • P7 P2 • ql Pi • q2 - 4 P5 • p, pi • ql p2 • q2 + 4 pz • Pv Pi • ql Ps • q2 + 
4 pi • P7 P2 • ql P5 • q2 + 4 P2 • Ps Pv • ql Pi • q2 + 4 pi • ps pv • ql pz • q2 + 
4 p2 • P7 Pi • ql P5 • ql + 4 pi • P7 P2 • ql P5 • ql + 4 P2 • P5 Pi • q2 p7 • q2 + 

4 Pi • P5 P2 • q2 P7 • q2 + Ps • P7 - 2 S Pl • P7 P2 • P5 - 2 S Pl • P5 P2 • P7 ) + 

1 6 P3 • P4 P6 • P8 P5 • ql P7 • q2 (P2 • ql Pi • q2 + pi • ql p2 • q2 ) + 
4 s tmass'' pa • P4 Ps • Pi Pe • Ps 

SetSF; 

AbsoluteTimingj^tmp2 = 

Expand FeynCalcFormLink[formexample, Print ^ False] /. |sP[ql, ql] ^tmass^, 



SP[q2, q2] -»tmass^, SP[pl, p2] -» — -emass^, SP[ql, q2] ^ — - tmass' 




{0. 167 021, 4 s tmass^ SP [pi, p7] SP [p2, p5] SP [p3, p4] SP [p6, p8] + 
4 s tmass^ SP [pi, p5] SP [p2, p7] SP [p3, p4] SP [p6, p8] + 
4 emass^ s tmass^ SP [p3, p4] SP [p5, p7] SP [p6, p8] - 
2 tmass^ SP [p3, p4] SP [p5, p7] SP [p6, p8] + 
4 s tmass'' SP [p3, p4] SP [p5, p7] SP [p6, p8] + 

8 tmass^ SP [pi, q2] SP [p2, ql] SP [p3, p4] SP [p5, p7] SP [p6, p8] + 
8 tmass^ SP [pi, ql] SP [p2, q2 ] SP [p3, p4] SP [p5, p7] SP [p6, p8] - 
8 tmass^ SP [pi, ql] SP [p2, p7] SP [p3, p4] SP [p5, ql] SP [p6, p8] - 
8 tmass^ SP [pi, p7] SP [p2, ql] SP [p3, p4] SP [p5, ql] SP [p6, p8] - 
8 tmass^ SP [pi, ql] SP [p2, p7] SP [p3, p4] SP [p5, q2 ] SP [p6, p8] - 
8 tmass^ SP [pi, p7] SP [p2, ql] SP [p3, p4] SP [p5, q2 ] SP [p6, p8] - 
8 tmass^ SP [pi, q2] SP [p2, p5] SP [p3, p4] SP [p6, p8] SP [p7, ql] - 
8 tmass^ SP [pi, p5] SP [p2, q2 ] SP [p3, p4] SP [p6, p8] SP [p7, ql] - 
8 emass^ tmass^ SP [p3, p4] SP [p5, ql] SP [p6, p8] SP [p7, ql] - 
8 emass^ tmass^ SP [p3, p4] SP [p5, q2 ] SP [p6, p8] SP [p7, ql] + 
4 s tmass^ SP [p3, p4] SP [p5, q2] SP [p6, p8] SP [p7, ql] - 
8 tmass^ SP [pi, q2] SP [p2, p5] SP [p3, p4] SP [p6, p8] SP [p7, q2 ] - 
8 tmass^ SP [pi, p5] SP [p2, q2 ] SP [p3, p4] SP [p6, p8] SP [p7, q2 ] + 
8 emass^ s SP [p3, p4] SP [p5, ql] SP [p6, p8] SP [p7, q2 ] - 
8 emass^ tmass^ SP [p3, p4] SP [p5, ql] SP [p6, p8] SP [p7, q2] + 
4 s tmass^ SP [p3, p4] SP [p5, ql] SP [p6, p8] SP [p7, q2 ] + 
16SP[pl, q2] SP[p2, ql] SP[p3, p4] SP[p5, ql] SP[p6, p8] SP[p7, q2] + 
16SP[pl, ql] SP[p2, q2] SP[p3, p4] SP[p5, ql] SP[p6, p8] SP[p7, q2] - 
8 emass^ tmass^ SP [p3, p4] SP[p5, q2 ] SP [p6, p8] SP[p7, q2]} 
TraditionalForm[tmp2 // FullSimplify] 
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(tmass^ (/75 • is (-2 emass^ + s - 2 tmass^) - 4 p2-ql /?i • q2 - 4 • ql p2- q2) + 2 • ql (2 emass^ 
P5-ql+ (2 emass^ - s) ps ■ q2 + 2 p2 ■ ps Pi ■ q2) +4 p2- pi pi-ql (ps ■ ql + P5 ■ q2)) 
2 p-i ■ q2 (p5 • ql {2 emass^ [s - tmass^) + 4 p\ ■ q\ p2- (\2 + s tmass^) - 2 emass^ tmass^ pi-q2 + 

Pi ■ q2 (4 p2-ql P5-qi-2 tmass^ p2 ■ P5)) + 
2tmass^ Pi ■ pi (2 p2 ■ ql {ps ■ q\ + ps ■ q2) - s p2- ps) + 
2 tmass^ pi ■ ps (2 p2 ■ q2 (pj ■ ql + pT q2) - s p2 ■ pi)) 
Expand [tmp2 - R3] 




B. Double Bremsstrahlung in r leptonic radiative decay 

This example is provided by Matteo Fael, it refers to a part of the calculation done 
Initialization 

UplUplbar = (*((1 + GA[5] . GS[n])/2) . (GS[pl] + M) ; *)(GS[pl] + M) ; 
Up2Up2bar = (GS[p2] + m) ; 



Bremsstrahlung 



Amplitudes without spinors 
' Gf 



Ma = 



*GA[a] . (1-GA[5]) 



(I* (GS[pl-kl-k2] +M) / (-2*SP[pl, kl] -2*SP[pl, k2] +2*SP[kl, k2])) 

I (GS[pl-kl] +M) 'I 
(-I*e*GA[p]) . I* . (-I*e*GA[^J]) 



■2 * SP [pi, kl] 

y"{\-y^)Gf i(y.(-ki-k2 + pi)+M) 



Mb = 



'_G^ 



-2ki ■ pi -2k2- pi + 2ki ■ k2 



*GA[a] . (1-GA[5]) 



.(-ieyP). 



( i(y(p,-ky)+M )\ 
2ki- pi 



(I* (GS[pl-kl-k2] +M) / (-2*SP[pl, kl] - 2*SP[pl, k2] + 2*SP[kl, k2])) . 

/ (GS[pl -k2] +M) 
(-I*e*GA[M]) . I* . (-I*e*GA[p]) 



■2 * SP [pi, k2] 

y".{l-y')Gf i(y.(-k,-k2 + pi)+M) 



Mc = (-1 *e*GA[p] ) . 



.(-ieyO- 

2ki- p[-2k2- pi+2ki-k2 

(GS[p2 + k2] +m) 



iiyiPi-k2) + M)] 
2k2- pi 



.(-ieyP) 



I * 



Gf 



*GA[a] . (1-GA[5]) 



2 * SP [p2, k2] 

(GS[pl - kl] +M) 



-2 * SP [pi, kl] 



(-I*e*GA[^J]) 
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i-ieyP). . \.i-iey^) 

2k2-p2 yfY V 2ki-pi 



( (GS[p2 + kl] +m) '\ 

Md= (-I*e*GA[iU]) . I* . 

I 2*SP[p2, kl] J 

I (GS[pl-k2] +M) \ 

I * . (-1 *e*GA[p] ) 

V -2*SPrpl, k21 ) 



*GA[a] . (1 - GA[5] ) 



i{y{h+p2) + m) y"{\-^)Gf ( i(y . (p, -k2) + M)\ 

i-iey"). . U-ieyP) 

2ki-p2 ^pl \ 2k2-pi ) 

[ (GS[p2 + k2] +m) A 

Me = (-I*e*GA[p]) . I* . (-I*e*GA[iu]) . 

V 2 *SP[p2, k2] ) 

(I* (GS[p2 + kl + k2] +m) / (2*SP[p2, kl] + 2*SP[p2, k2] + 2*SP[kl, k2])) 
( Gf 



*GA[a] . (1 -GA[5] ) 

IVi" 

i{y-(k2 + P2)+m) i(y-(ki+k2 + p2) + m) y''.{l-y^)Gf 

i-ieyP). .(-teyO. . 

lk2- P2 2ki ■ p2 + 2k2- p2 + '2ki ■ k2 ^|2 

I (GS[p2 + kl] +m) ^ 

Mf = (-I*e*GA[Ai]) . I* . (-I*e*GA[p]) . 

V 2*SP[p2, kl] ) 

(I* (GS[p2 + kl + k2] +m) / (2*SP[p2, kl] +2*SP[p2, k2] + 2*SP[kl, k2])) 

Gf ^ 
— — *GA[a] . (1-GA[5]) 

i iy ■ {kl + p2) +m) i{y (ki +k2 + P2) + m) y"{l - Gf 

{-iey^). .{-ieyP). . 

2ki- p2 2ki- p2 + 2k2- P2 + 2ki-k2 ^JY 

Mbrem = Ma + Mb + Mc + Md + Me + Mf ; 

Mbremstar = FCE [ComplexCon jugate [Mbrem] / . {a -» /3, ix -* v, p -» ct} ] ; 
ComplexCon jugate [Ma] 

Gf r^Xr • iPi - h) +M).yP.{y ■ {-k, - ^2 + /^i) + M){y' + l).y" 

2yp2 k\ ■ p\ (-2ki ■ p\-2k2- Pi+2k\-k2) 
Squared Amplitude 

ATTENTION! The tensor of Neutrini changes in the double photon emission! 
1 

Neutriniyy = * ((FV[pl, a] -FV[p2, a] -FV[kl, a] -FV[k2, a]) * 

3 *Pi 

(FV[pl, li] -FV[p2, ^] -FV[kl, li] -FV[k2, ^] ) - 
MT[a, J3] * (M'^2 +m'^2 - 2 * SP [pi, p2] -2*SP[pl, kl] + 

2*SP[p2, kl] - 2*SP[pl, k2] +2*SP[p2, k2] + 2*SP[kl, k2])) 

—((-^1" - ^2" + Pf - Pi") - A:2^ + Pi^ - P2^) - 

g"'^[-2ki ■ pi+2ki-p2-2k2- pi+2k2- p2 + 2ki-k2+m^ -2pi- P2)) 
DlracSimplify[DiracGainma[7] , DiracSubstltute67 -* True] 
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2 2 

Fbrem = DiracTrace[MT [ju, v]MT[p, a] Mbrem.UplUplbar.Mbreinstar.Up2Up2bar] /. 
{ (1 + GA[5] ) ^ (2 GA[6] ) , (1 - GA[5] ) ^ (2 GA[7] ) , (1 - GA[6] ) ^ (GA[7] ) } 



tr 



i (m + y • (fei + p2)) i{m + y-(kx+k2 + pi)) r"-(2 J^) Gf 



Iky p2 2ki-k2 + 2ki- p2 + 2k2- P2 yjY 

i{m + y(ky+ p2)) y"{2 y'^]Gf ( i{M+y(px- ^2)) 



lki-p2 



2k2-pi 



{-ieyP) + 



iim + y-{k2 + P2)) i{m + y-iki+k2 + P2)) 7"-{^y^)Gf 
i-ieyf). i-iey^). . + 



(-ieyP). 



2k2- P2 2k\-k2 + 2k\- P2 + 2k2- P2 ypl 

iim + y ■ik2 + P2)) r"-(2 7')Gf ( i{M+y (p, - k,)) 



2k2-p2 

r{2y')Gf i(M+yi-ki-k2 + pi)) 



Y2~ 2ki-k2-2ki- pi-2k2- pi 

ri2y')Gf i(M +yi-ki-k2 + pi)) 
yj~2 2k\ ■k2 -2k\ • p\ -2k2- p\ 



.{-ieyi") 
.(-ieyf). 



2ki-pi 

iiM +yipi-k2)) 



2k2-p\ 

( i(M + y-(pi-ki)) 
2k\ ■ pi 



i-iey^) + 

(-ieyP) + 



.{-ieyi") 



{M +y-px). 



-{y\{M +y[pi- h)).y''m + y ■ {-ki -k2 + pi)).{2 y^).yP Gf e") j 



{2yf2ki-pii2ki-k2-2ki-pi-2k2-pi))- 

rm + r ■ (/'I - fe)).(2 /)./.(m + r • + Vi))-Y Gf 

4^fY ki- p2k2- pi 

(r".(M +y(pi- k2)).Y.{M +y(-h-k2 + p,)).[2y%yP Gf e") 

{2 yp2 (2ki ■ k2-2ki ■ pi -2k2- pi) k2 ■ Pij - 

Y.{M +yipi- k,)).{2 /)y .(m + y ■ + Pi))-?^ Gf 
+ 

4^fYkl■plk2■p2 
((2 /)./.(m + yiki+k2 + P2)).y^.im + yiki+ p2)).Y Gf e^) / 

{2yf2 kyp2{2ki-k2+2kyp2 + 2k2- P2)] + 
((2 /).7^.(w + y{ki+k2 + P2)).f.{m + y{k2 + P2)).y'' Gf e") 

{2y[Y k2- P2i2ki-k2 + 2ki-p2 + 2k2- P2)) 



\ \ 
.im + yp2)g'"gf"' 



AbsoluteTiming[Amp = FeynCalcFormLink[Neutrini°2CY * Fbrem, 
IDStatements -* 

"id kl.kl = 0; id k2.k2 = 0; id pl.pl = M'^2; id p2.p2 = m^2;" 

];] 
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Symbols e,Gsubf,m,M; 
Vectors kl,k2,pl,p2; 
AutoDeclare Index lor; 
Format Mathematica; 

L resFL = ((d_(lor3,lor4)*d_(lor5,lor6)*(-((e^2*Gsubf*g_(l,lorl)*(g7_(l)/2)*(M-g_(l,kl) 
-g_(l,k2)+g_(l,pl))*g_(l,lor5)*(M-g_(l,kl)+g_(l,pl))*g_(l,lor3))/(sqrt_(2) 
*pl(kl)*(2*k2(kl)-2*pl(kl)-2*pl(k2))))-(e^2*Gsubf*g_(l,lorl)*(g7_(l)/2) 
*(M-g_(l,kl)-g_(l,k2)+g_(l,pl))*g_(l,lor3)*(M-g_(l,k2)+g_(l,pl)) 
*g_(l,lor5))/(sqrt_(2)*(2*k2(kl)-2*pl(kl)-2*pl(k2))*pl(k2))-(e^2*Gsubf 
*g_(l,lor3)*(m+g_(l,kl)+g_(l,p2))*g_(l,lorl)*(g7_(l)/2)*(M-g_(l,k2) 
+g_(l,pl))*g_(l,lor5))/(2*sqrt_(2)*pl(k2)*p2(kl))-(e^2*Gsubf*g_(l,lor5) 
*(m+g_(l,k2)+g_(l,p2))*g_(l,lorl)*(g7_(l)/2)*(M-g_(l,kl)+g_(l,pl)) 
*g_(l,lor3))/(2*sqrt_(2)*pl(kl)*p2(k2))+(e^2*Gsubf*g_(l,lor3)*(m+g_(l,kl) 
+g_(l,p2))*g_(l,lor5)*(m+g_(l,kl)+g_(l,k2)+g_(l,p2))*g_(l,lorl) 
*(g7_(l)/2))/(sqrt_(2)*p2(kl)*(2*k2(kl)+2*p2(kl)+2*p2(k2)))+(e^2*Gsubf 
*g_(l,lor5)*(m+g_(l,k2)+g_(l,p2))*g_(l,lor3)*(m+g_(l,kl)+g_(l,k2)+g_(l,p2)) 
*g_(l,lorl)*(g7_(l)/2))/(sqrt_(2)*p2(k2)*(2*k2(kl)+2*p2(kl)+2*p2(k2)))) 
*(M*gi_(l)+g_(l,pl))*(-((e^2*Gsubf*g_(l,lor4)*(M-g_(l,kl)+g_(l,pl)) 
*g_(l,lor6)*(M-g_(l,kl)-g_(l,k2)+g_(l,pl))*(g6_(l)/2)*g_(l,lor2))/(sqrt_(2) 
*pl(kl)*(2*k2(kl)-2*pl(kl)-2*pl(k2))))-(e^2*Gsubf*g_(l,lor6)*(M-g_(l,k2) 
+g_(l,pl))*g_(l,lor4)*(M-g_(l,kl)-g_(l,k2)+g_(l,pl))*(g6_(l)/2) 
*g_(l,lor2))/(sqrt_(2)*(2*k2(kl)-2*pl(kl)-2*pl(k2))*pl(k2))-(e^2*Gsubf*g_(l,lor6) 
*(M-g_(l,k2)+g_(l,pl))*(g6_(l)/2)*g_(l,lor2)*(m+g_(l,kl)+g_(l,p2)) 
*g_(l,lor4))/(2*sqrt_(2)*pl(k2)*p2(kl))-(e^2*Gsubf*g_(l,lor4)*(M-g_(l,kl) 
+g_(l,pl))*(g6_(l)/2)*g_(l,lor2)*(m+g_(l,k2)+g_(l,p2))*g_(l,lor6))/(2 
*sqrt_(2)*pl(kl)*p2(k2))+(e^2*Gsubf*(g6_(l)/2)*g_(l,lor2) 
*(m+g_(l,kl)+g_(l,k2)+g_(l,p2))*g_(l,lor6)*(m+g_(l,kl)+g_(l,p2)) 
*g_(l,lor4))/(sqrt_(2)*p2(kl)*(2*k2(kl)+2*p2(kl)+2*p2(k2)))+(e^2*Gsubf 
*(g6_(l)/2)*g_(l,lor2)*(m+g_(l,kl)+g_(l,k2)+g_(l,p2))*g_(l,lor4) 
*(m+g_(l,k2)+g_(l,p2))*g_(l,lor6))/(sqrt_(2)*p2(k2)*(2*k2(kl)+2*p2(kl) 
+2*p2(k2))))*(m*gi_(l)+g_(l,p2))*((-kl(lorl)-k2(lorl)+pl(lorl)-p2(lorl))*(-kl(lor2) 
-k2(lor2)+pl(lor2)-p2(lor2))-d_(lorl,lor2)*(m^2+M^2+2*k2(kl) 
-2*pl(kl)-2*pl(k2)+2*p2(kl)+2*p2(k2)-2*p2(pl))))/(3*pi_)); 

trace4,l; 

contract 0; 

.sort; 

id kl.kl = 0; 
id k2.k2 = 0; 
id pl.pl = M^2; 
id p2.p2 = m^2; 
.sort; 

#call put("%E", resFL) 
#fromexternal 

Piping the script to FORM and running FORM 

Time needed by FORM : 4.583 seconds. FORM finished. Got the resuh back to Mathematica as a string. 
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Start translation to Mathematica / FeynCalc syntax 

Total wall clock time used: 13.28 seconds. Translation to Mathematica and FeynCalc finished. 
{13.378199, Null} 
Variables [Amp] 



% 1 1 InputForm 

{e, m, M, SP[kl, k2], SP[kl, pi], SP[kl, p2], 

SP[k2, pi], SP[k2, p2], SP[pl, p2], Subscript[G, f ] } 
FORMAmp = Expand [Amp] ; 
Length @ FORMAmp 
2612 

Check with FeynCalc 

ScalarProduct[kl, kl] =0; 
ScalarProduct[k2, k2] =0; 
ScalarProduct[pl, pi] = M^2; 
ScalarProduct[p2, p2] =m^2; 

The naive approach is too slow. 

(*AbsoluteTiming[FCAmp = TR [Neutriniyy * Fbrem];]*) 

Do this instead (also still quite slow, but good enough for checking equality) 
AbsoluteTiming FCAmpl = Contract [Neutriniyv DiracGammaExpand[FCI [Fbrem] ] ] / . 



DiracTrace ^ DiracTrick; DotExpand2 [expr_] := 

expr/. Dot -» Hold [Dot] //. {Hold[Dot] [a , b_Plus, c ] m 

(Distribute[Hold[Dot] [a, b, c] ] //. {Hold[Dot] [aa , bb_ cc_, dd ] 



{Hold[Dot] [] :^1, Hold[Dot] :^Dot}; Print [ 
"noncommutative expansion, time needed: ", 
AbsoluteTiming [ 
FCAmp2 = Expand[DotExpand2 [FCAmpl] , Dot] ; ] ] ; 
Print [ "collecting traces , time needed: ", 
AbsoluteTiming [ 

FCAmp3 = Collect2 [DiracTrick[FCAmp2] , Dot, Factoring-* False] ; ] ] ; 
Print ["Length [FCAmp3] = ", Length [FCAmpS] ] ; 




bb Hold [Dot] [aa, cc, dd] /; NonCommFreeQ[bb] , Hold [Dot ] [aa 



bb_, dd ] :-»bb Hold [Dot] [aa, dd] /; NonCommFreeQ[bb] } ) } /. 



Print "time doing the traces 



Length [ FCAmpS ] 




AbsoluteTiming 



zeit = AbsoluteTime [ ] ; FCAmp4 = 




0. 



Print Temporary [ { i , 



, Round [AbsoluteTime [] -zeit]}] 



Expand [ExpandScalarProduct@TR [FCAmpS HiJ ] ] 




FCAmpS = Expand[ScalarProductExpand[FCAmp4] ] ; 
Print ["length of FCAmpS = ", Length [FCAmpS] ] 

noncommutative expansion, time needed: {44.751183, Null} 
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collecting traces , time needed: {169.767466, Null) 

Length[FCAmp3] = 694 

time doing the traces {62.221846, Null} 

length of FCAmp5 = 1896 

{277.372532, Null} 

AbsoluteTiming[dif f = Collect [FCE [FCAmpS - Amp] , M, Factor]] 

{16.162748, 0) 

Tlie results of FORM and FeynCalc agree, but FORM is much faster. There is a some 
noticeable overhead in translating the returned string to FeynCalc. In a future version 
of FormLink we will improve this by either using parallel features of Mathematica or by 
writing a special MathLink program similar to FormGet . tm. 

VI. Summary 

We have implemented FormLink and FeynCalcFormLink, two MathLink-hased pro- 
grams implemented in Mathematica, C and FORM, which embed FORM in Mathematica 
and FeynCalc. A non-MathLink-hased file-based function RunForm has been also put 
into the FormLink package. For both packages we provide configurable functions for the 
syntax conversions. A limited subset of the syntaxes from both programs are automatically 
translated to each other by our two Mathematica packages. In this way we can combine the 
speed of FORM with the generality of Mathematica. A simple installation facility has been 
provided. 

We hope to be able to extend the functionality and range of applicability of the packages 
in the future. 

VII. Licenses 

FormLink and FeynCalcFormLink are covered by the GNU Lesser General Pubhc 
License, like FORM. The conditions for the use of FORM are laid out here: |http: // 
www_. nikhef .nl7~f orm/license/l icense .html and should be followed of course also when 
using FormLink/FeynCalc FormLink. The license for using Mathematica is given here: 
|http : / / www . wolfram . com/legal/ agreements/wolf ram-mathematica . html, 
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